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Abstract 

The short-time asymptotic behavior of option prices for a variety of models with jumps has received much 
attention in recent years. In the present work, a novel second-order approximation for ATM option prices under the 
CGMY Levy model is derived and, then, extended to a model with an additional independent Brownian component. 
Our method of proof is based on an integral representation of the option price involving the tail probability of the 
log-return process under the share measure and a suitable change of probability measure under which the process 
becomes stable. This approach is sufficiently efficient to produce the third-order asymptotic behavior of the option 
prices and, moreover, is expected to apply to many other popular classes of Levy processes which satisfy the 
fundamental property of being stable under a suitable change of probability measure. Our results shed new light on 
the connection between both the volatility of the continuous component and the jump parameters and the behavior 
of ATM option prices near expiration. In the case of an additional Brownian component, the second-order term, 
in time-f, is of the form d2 t^^~^^^^, with the coefficient d-z depending only on the overall jump intensity parameter 
C and the tail- heaviness parameter Y. This extends the known result that the leading term is {a/^/2Tv)t''^^ , where 
a is the volatility of the continuous component. In contrast, under a pure-jump CGMY model, the dependence on 
the two parameters C and Y is already reflected in the leading term, which is of the form dit^^^ . Information on 
the relative frequency of negative and positive jumps appears only in the second-order term, which is shown to be 
of the form d2t and whose order of decay turns out to be independent of Y. The third-order asymptotic behavior 
of the option prices as well as the asymptotic behavior of the corresponding Black-Scholes implied volatilities are 
also addressed. Our numerical results show that first-order term typically exhibits extremely poor performance and 
that the second-order term significantly improves the approximation's accuracy. 
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1 Introduction 

It is generally recognized that the standard option pricing model of Black-Scholes is inconsistent with options data, while 
remaining a widely used model in practice because of its simplicity. Exponential Levy models generalize the classical 
Black-Scholcs setup by allowing jumps in stock prices while preserving the independence and stationarity of returns. 
There are several reasons for introducing jumps in financial modeling. First of all, asset prices do jump, and some 
risks simply cannot be handled within continuous-paths models. Second, historical asset prices exhibit distributions 
with so-called stylized features, such as heavy tails, high kurtosis, volatility clustering and leverage effects, which are 
hard to replicate within purely-continuous frameworks. Finally, market prices of vanilla options exhibit skewed implied 
volatilities (relative to changes in the strikes), in contrast to the classical Black-Scholes model which predicts a flat 
implied volatility smile. Moreover, the fact that the implied volatiliy smile and skewness phenomenon becomes much 
more pronounced for short maturities is a clear indication of the presence of jumps. 
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One of the first apphcations of jump processes in financial modehng is due to Mandelbrot [SBj . who suggested a pure- 
jump stable Levy process Z to model power-like tails and self-similar behavior in cotton price returns. Merton [28j and 
Press [31] subsequently considered option pricing and hedging problems under an exponential compound Poisson process 
with Gaussian jumps and an additive independent non-zero Brownian component. A similar exponential compound 
Poisson jump-diffusion model was more recently studied in Kou (22j . where the jump sizes are distributed according to 
an asymmetric Laplace law. For infinite activity exponential Levy models, Barndorff-Nielsen introduced the normal 
inverse Gaussian (NIG) model, while the extension to the generalized hyperbolic class was studied by Eberlein, Keller 
and Prause [10]. Madan and Seneta [25] introduced the symmetric variance gamma (VG) model while its asymmetric 
extension was later studied by Madan and Milne [21] and Madan, Carr and Chang [23]. Both models are built on 
Brownian subordination; the main difference being that the log-return process in the NIG model is an infinite variation 
process with stable like (a = 1) behavior of small jumps, while in the VG model, the log-price is of finite variation 
with infinite but relatively low activity of small jumps. The class of "tempered stable" processes was first introduced 
by Koponen [21] and further developed by Carr, Geman, Madan and Yor [5] , who introduced the terminology CGMY. 
The CGMY model is a particular case of the more general KoBoL class of [1] and was also previously proposed for 
financial modeling in [7J and [27]. Nowadays, the CGMY model is considered to be a prototype of the general class of 
models with jumps and enjoys widespread applicability. 

Stemming in part from its importance for model calibration and testing, small-time asymptotics of option prices 
have received a lot of attention in recent years (see, e.g., [2], [3], [U], [12], [13], [16], [II], [^, [19], [20], [30], [32], 
[55]). We shall review here only the studies most closely related to ours, focusing in particular on the at-the- money 
(ATM) case. Carr and Wu [S] first analyzed, partially via heuristic arguments, the first order asymptotic behavior of 
an Ito semimartingale with jumps. Concretely, ATM option prices of pure-jump models of bounded variation decrease 
at the rate 0(t), while they are just 0{\/t) under the presence of a Brownian component. By considering a stable 
pure-jump component, [9] also showed that, in general, the rate could be 0(i'^), for some /3 e (0,1). Muhle-Karbe 
and Nutz |29j formally showed that, under the presence of a continuous-time component, the leading term of ATM 
option prices is of order y^, for a relatively general class of Ito models, while for a more general type of Ito processes 
with a-stable-like small jumps, the leading term is 0{i}l°') (see also [131 Proposition 4.2], [TJl Theorem 3.7], and [3H1 
Proposition 5] for related results in exponential Levy models). However, none of the these papers obtained second or 
higher order asymptotics for the ATM option prices, which are arguably more relevant for calibration purposes, given 
that the most liquid options are of this type. 

In the present paper, we study the small-time behavior for at-the-money (ATM) call (or equivalently, put) option 
prices 

E(5t-5o)+ = 5oE(e^' -1)+, (1.1) 

under the exponential Levy model 

St-.^S^e''^ (1.2) 
where X is the superposition of a CGMY Levy process (Lt)t>S) ^nd of an independent Brownian motion (crWt)t>o; i-e., 

Xt = U^cjWt. (1.3) 

where (M^t)t>o is a standard Brownian motion independent of L. Here, as usual, is the positive part of x. The first 



order asymptotic behavior of (1.1) in short-time under the model ( |1.3[ ) takes the form: 

\\Tiit-^l^mSt - SqY = SoE{Z+), (1.4) 

where Z is a symmetric stable random variable with a = Y under P. When a 0, Z ^ A/'(0, cr^) {a — 2) and, thus, 
E(Z+) = a/^/2Tr (see [3S] and [32])- When cr = and a — Y, the characteristic function of Z is explicitly given (see 
[T3] and (3S1) by 

Ee'"^ ^ g-2cr(-y)|cos(iy7r)||«r^ 



In that case, (see (25.6) in [37]), 

di := E{Z+) = (^1 - (^2CT{-Y) 



cos 



Try 



l/Y 

(1.5) 
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Interestingly enough, under the presence of a continuous component, the first-order asymptotic term only reflects 
information on the continuous-time volatility, in sharp contrast with the pure-jump case where the leading term 
depends on the overall jumps-intensity parameter C and the index y, which in turn controls the tail-heaviness of the 
distributions. 



Below, we also obtain a second order correction term for the approximation (1.4). The derivation of the second- 
order results builds on two facts. First, as in [T^, we make use of the following representation of Carr and Madan 

—E{St-So)+=¥*iXt>E)= e-''¥*iXt>x)dx, (1.6) 
•^0 Jo 

where P* is the martingale probability measure obtained when one takes the stock as the numeraire (i.e., V*{A) := 
E(5'(1a)) and E is an independent mean-one exponential random variable under P*. The measure P* is sometimes 



called the share measure (see [6]). Notice that under P*, {Xt)t>o also admits a decomposition similar to (1.31 



Xt^L*+aWt*, i>0, (1.7) 

where W* := {Wf)t>o is a Wiener process and L* :— (L()t>o is also a CGMY process, independent of W* . Second, 
we change probability measures from P* to a probability measure P, under which (L^ )t>o is a stable Levy process and 
{Wf)t>o is still a standard Brownian motion independent of L* . We show that the second-order asymptotic behavior 



of the ATM call option price (1.1) in short-time is then of the form 

^E{St- So)+ =dit^ +d2t + oit), {t^O), 
Jo 

in the pure-jump CGMY case (cr = 0), while in the case of a non-zero independent Brownian component (cr 7^ 0), 

^^{St - So)+ = dit^ + d2t^ + o (t^) , {t ^ 0), 

for different constants di and ^2 that we will determine explicitly. To wit, we found that, under the presence of a 
nonzero Gaussian component, the second-order term depends only on the overall jump intensity parameter C and the 
tail-heaviness parameter Y . The parameters G and M (which control the relative frequency of negative and positive 
jumps) do not appear until the next order term. However, for a pure-jump case, the parameters G and M are already 
present in the second-order term. The above asymptotic behaviors should also be compared to the corresponding 
behavior under the standard Black-Scholes model, where it is known that (see, e.g., [TBI Corollary 3.4]) 

E(e-^*-^* _ 1)+ = _ -^^t^- + 0[ti). 

Our method of proof is sharp enough to produce the third-order asymptotic behavior of the option prices (see Remark 



3.4 and 4.4 below) and, moreover, is expected to apply to other popular classes of Levy processes, which satisfy the 
fundamental property of being stable under a suitable change of probability measure such as tempered stable processes 
in the sense of Rosihski [35] (this will be presented elsewhere). Finally, the asymptotic behavior of the corresponding 
Black-Scholes implied volatilities are also addressed. 

The present paper is organized as follows. Section 2 contains preliminary results on the CGMY model, some 
probability measure transformations, and asymptotic results for stable Levy processes which will be needed throughout 
the paper. Section 3 establishes the second-order asymptotics of the call option price under the pure-jump CGMY 
model (cr — 0). Section 4 establishes the second-order asymptotics of the call-option price under the CGMY model 
with an additional independent non-zero Brownian component (cr ^ 0). In Section [sj we assess the performance of 
our asymptotic expansions through a detailed numerical analysis. The proofs of our main results are deferred to the 
Appendices. 
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2 Setup and preliminary results 



2.1 The CGMY model 

Throughout, {Lt)t>o stands for a CGMY Levy process defined on a complete filtered probability space (il, (J^t)(>o, P) 
with corresponding parameters C, G, M > and Y G (I, 2). That is, i is a pure-jump Levy process with characteristic 
function 

E (e™-^*) = exp (t [icu + CT{-Y) ((M - iuf + (G + luf - - G^)]) . (2.1) 

Let Xt=i7Wt + Lt, t > Q, where {Wt)t>o is a standard Brownian motion, independent of {Lt)t>Oj defined on 
(fi, J", (J"t)t>o,P). We call the process X the (generalized) CGMY model. 

We assume zero interest rate and that P is a martingale measure for the exponential Levy model 5*4 = S'oe"^'. In 
particular, M > 1 and the characteristic function ipt of Xf is given by 



(pt{u) =E(e™^') =exp (i 

with 



+ Gr(-r) ((Af - iuf + (G + tuf - - G^) 



(2.2) 



2 

c = -Gr(-F) ((M - I)^ + (G + 1)^ - - G^) - y ; (2.3) 

see, e.g.. Proposition 4.2 in ^38^. In particular, we note that 7 := EXi = ELi is given by 

7 = c-Grr(-y)(M^-i-G^-i). (2.4) 

The Levy triplet of {Xt)t>o (relative to the truncation function 2:l{|j.|<x}) is denoted by (6, (t^,i/). Thus, v and b are 
given by 

^(^^) = y ^i+Y l{->0} + ^^TTF l{-<0} j '^a:, (2.5) 

6 = c- / a;i/(da;)-Gyr(-y)(Af^-i-G^"^). (2.6) 
J|x|>i 

Without loss of generality, we also assume throughout that {Xt)t>o is the canonical process Xt{uj) — uj{t) defined 
on the canonical space fl = D([0, oo),IR) (the space of cadlag functions lo : [0, oo) — >■ K) equipped with the cr-field 
J- = a{Xs : s > 0) and the right-continuous filtration Tt '■— HgytO'iXu : u < s). 



2.2 Probability measure transformations 

Following a density transformation construction as given in Sato |37j (see Definition 33.4 and Example 33.4 there) and 
using the martingale condition Ee"^' = 1, we define a probability measure P* on (il, J^) such that 



{t > 0); 



(2.7) 



i.e., ¥*{B) — E (e^^'ls), for any B ^ Ft and t > Q. The measure P* can be interpreted as the martingale measure 
when using the stock price as the numeraire. Under P*, {Xt)t>o is also a Levy process and its characteristic function 
is given by 



E*(e^"-^') =exp ( t 
with (see Appendix [C|) 



+ Gr(-y) ((Af* - lu)^ + {G* + luY - M*^ - G*^'^ j 



(2.8) 



Af*=Af-l, G*=G + 1, c=c + (j^ 
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It is clear from (2.8) that, under P*, {Xt)t>o can also be decomposed as in (1.7), where {Wf)t>o is again a Wiener 
process while {Ll)t>o is still a CGMY process, independent of W* , but with parameters C, Y, M — M* , and G = G* . 
Hereafter, we denote the Levy triplet of {Xt)t>o under P* by (6*, (cr*)^,z/*), where a* = a, i'*(dx) = e^v{dx)^ and 

b* := c* - / XV* [dx) - CYT{-Y){{M*Y-^ - (G*)^-i). (2.9) 

As explained in the introduction, an important tool in the sequel is to change the probability measures from P* to a 
probability measure P, under which (Lj )f>o is a stable Levy process and (VFt*)t>o is still a Wiener process independent 
of L*. Concretely, let 

v{dx) ■.^C\x\^^-^dx, b^b*+[ x{D-v*){dx). 

J\x\<l 

Note that v is the Levy measure of a symmetric stable Levy process and, also, 

v{dx) = e^(^V*(dx), 

with 

tf{x) := M*x l{x>a} - G*x l{x<o}- 

Hence, by virtue of Theorem 33.1 in [37], there exists a probability measure P locally equivalent]^ to P* such that 
(Xt)(>o is a Levy process with Levy triplet (6, cr^, v) under P. Throughout, E denotes the expectation under P. 

In light of (2.9) and since EXi = EL| = b + S{\x\>i} xi'idx), it can be shown (see Appendix |c| that 



7 := EXi = -Cr{-Y) {{M - if + {G + if - - G^) + (2.10) 

Next, we recall that the centered process {Zt)t>o, given by 

Zt:=L;-tj, (2.11) 
is symmetric and strictly F-stabl^ under P and, thus, is self-similar; i.e., 

(t Zut)u>o = iZu)u>o, (2-12) 
for any t > 0. We also need the following representation of the density process (see Theorem 33.2 in |37)): 



dP- 

with 

Uf :~ lim 



e^S (2.13) 



( ^ ^{AX^)-tf {e^^^^-l)v*{dx)] 

\s<t:\AXs\>e "'l^^^ / 



The process {Ut)t>o can be expressed in terms of the jump-measure N{dt,dx) :— #{(s, AX^) E dt x dx} of the 
process {Xt)t>o and its compensator N{dt,dx) := N{dt,dx) — v{dx)dt (under P); namely, 

Ut = M*C/t+ - G*C7r + ?7t, (2.14) 



'^Equivalently, there exists a process {Ut)t such that P(i?) = E*(e^*lB), for t > and B £ J-t. 
■^Concretely, its scale, skewness, and location parameters are 2CT(—Y)\ cos(7ry/2)|, 0, and 0, respectively. 
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where 



UfT := I / xN{ds,dx), 

Iq "'(0,00) 



xN{ds, dx), 



"'(-00,0) 



ry.^C {e 



-M*x 



l + M*x)x-''-'dx + C / {e^''' - l-G*x)\x\-''-'dx 



= CT{-Y) {{M*f + {G*f) 



(2.15) 



(2.16) 



where in the last equaUty we used the analytic continuation presented in [37] (see (14.19) therein). Finally, let us also 
note the following decomposition of the process X in terms of the previously defined processes: 



Xt = Zt+t^ + (jW* = C/+ + U;- + i7 + aW*. 
The following table summarizes the different probability measures used in this paper: 



Prob. Measure 


{Lt)t distribution 


Density wrt P 


P 


CGMY(C, G, M, Y) 


1 


P* 


CGMY(C, G*,M*,F) 




P 


Stable(/3 = 0; a = V) 





(2.17) 



2.3 Some needed properties of stable Levy processes 



Let us now collect some well-known results on stable Levy processes needed in the sequel. First, from ^.15[ ), it is 
clear that {Uf^)t>o and {—U^)t>o are independent and identically distributed one-sided y-stable processed^ under P. 
Hence, the common transition density of Uj^ and —Uj~, denoted by p{t, x), exists (cf. [371 Proposition 2.5]). Moreover, 
the following result for the asymptotic behavior of the transition density is known (see, e.g., f35 and [H]): 

lmijp{t,u)^s{u), (UT^O), 

where s is the Levy density of the Levy process {U^)t>o. In particular, since by construction the Levy measure of 
{U^)t>o is I'+idu) = Cu~^~^l{u>o}du, the Levy density s is just Cw~^~^l{„>o}, u 7^ 0, and we get: 



t " 1{«>0} , u^O. 



(2.18) 



Equivalently, by the self-similarity of {Uf^)t>o, we have 



and (2.181 can be casted as follows by setting x = t u: 

p{\,x) - Cx'^~^, 



(x — >■ QO). 



As a consequence. 



C_ 



(2.19) 



(2.20) 



^Concretely, its scale, skewness, and location parameters are C| cos(7ry/2)|r(— y), 1, and 0, respectively. 
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as a; — > oo. Equivalently, plugging x = t v and using the self-similarity of and U , we recover the well-known 
result: 



lim ^Pf 

t->-o t \ 



In particular, there exists iV > 0, such that for all < t < 1 and v > satisfying t v > A^, 



2C 



p(;7+ > t-^^^v) < ^tv' 



-l/Y 



, 2C Y 



(-C/f > t 

The following result sharpens (2.22 1. Its proof is presented in Appendix [C| 

Lemma 2.1. There exists a constant < k < oo such that for any < t < 1 and u > 0, 

P(C7+ > t-^'^v) < KtV-^ , H-Ul > t-^'^v) < Ktv'^. 

Therefore, since Zi — — f/j", 

P(Zi > t-^/^v) = FiZt >v)< 2^+^Ktv-^ < Sutv-^, 
for any < t < 1 and v > 0. Note also that 

P{Zt >v)^ ¥{Zi > t-^/^v) - tv{[v, oo)) = t yv-^, (t 0), 

and that the probability density pz of Zi is such that 

Pz{v) Cv-^-'^, (w^oo). 
The following identity for Ut '■— M*U^ — G*U^ will also be needed in sequel: 



E 



(e-*^^"^^) =E*(e-*^^"^^*P.")r(e*^'"«'^^-) =exp(,yt). 



(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



The relation above follows from the representations (2.151, the independence of t/+ and U , and the form of the 
characteristic function of a Poisson integral. 



3 The pure-jump CGMY model 



In this section, we find the second-order asymptotic behavior for the at-the-money call option prices (1.1) in the pure 



jump CGMY model. The proofs of all results in the section are deferred to the Appendix]^ Throughout this section, 
{Xt)t>Q is a Levy process with triplet (6,0,:/) as introduced in Section [2j As explained in the introduction, the first 
order asymptotic behavior is given by (1.4). Before stating our first result, we need to rewrite the call option price 



(1.1 ) in a suitable form. 



Lemma 3.1. In terms of the probability measure P defined in (2.13) and the parameter 7 defined in (2.10), we have 
that 



1 1 



JO 



(e * ' ^^l{Zi>v}^ dv, 



(3.1) 



where Ut '.= M*U^ — 0*11^ , and {Ut)t>o and {U^ )t>Q are defined as in ('2.15) 



The following result gives the second-order asymptotic behavior of at-the-money call option prices under the pure- 
jump CGMY model. 
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Theorem 3.2. Under the exponential CG MY model (1.2) without Brownian component, 



lim t Y ' 



1 



where rj and 7 are respectively given as in (21^ and (21^ , and 

d -CT(-y)((M* + 1)^ + (G*)^). 



(3.2) 



(3.3) 



Remark 3.3. Using (3.3), (2.161, and (2.101, it turns out that 



^2 • = 



+ ^^LXl ((M - 1)^ - -{G+ 1)^ + G^) 



Hence, the second-order asymptotic behavior of the ATM cah option price (1.1) in short-time is 

1 



^0 



E{St~ Sa)+ =dit^ +d2t + oit), [t^O], 



(3.4) 



with di — E(^i ) given as in (1.5 1: 



^ ir (1 - 1 ) ( 2cr(-r) 



ttY 



Broadly speaking, the first-order term synthesizes only the information on the tail-heaviness index Y and the overall 
jump-intensity parameter G, while the second-order term incorporates also the information on the relative intensities 
of negative and positive jumps (controlled by the parameters G and M). Note also that ^2 < — Gr(— F) < — 2G. 



Remark 3.4. The proof of Theorem |3.2| (see Appendix |A|) also provides the higher order asymptotics of the ATM 
call option price under the pure jump CGMY model. Indeed, it is clear that the term D2 defined in ( A.5[ ) is such 
that D2{t) ^ —riE{Z^)t. Moreover, the second-order term of Di{t) therein can be shown to be 0{t^~^), while 
1)3 (t) ^ o(D2(i)), ast^O. Therefore, as t ^ 0, and since 1 < F < 2, 



So 



E{St - S'o)+ = dit^ + d2t - riK{Z+)t 



o{t 



(3.5) 



Let a{t) denote the ATM Black-Scholes implied volatility at maturity t with zero interest rates. The following 
result gives the asymptotic behavior of a-{t) as t — > 0. 



Proposition 3.5. Under the exponential CGMY model [1.2) without Brownian component, the implied volatility a{t) 
has the following small-time behavior: 



a{t)=aitY 2+0-2^2+0(^2), 



t ^ 0, 



where 



0-1 
0-2 



2ttE(Z{), 

y|Gr(-y) ((Af - 1)^ - - (G + 1)^ + G^) 



(3.6) 

(3.7) 
(3.8) 



4 The CGMY model with Brownian component 



In this part, we consider the CGMY model with non-zero Brownian component. Concretely, throughout, (Xt)tyf) is a 



Levy process with triplet (6, cr^, v) as introduced in Section [2] and a ^ 0. In that case, it follows from (2.2) that 



limE* (^exp{iuXt / Vi)^ = exp 



- 2 2 
-an 



8 



and, thus, [Xt/^/t) converges weakly to the centered Gaussian distribution with variance ct'^. Equivalently, recalhng 
that under P*, (H/j*)t>o is a standard Brownian motion, it fohows that 

(4.1) 



HmP*(Xf/Vt > x) = P*(crW^i* > x). 

t-i-O 

The first order asymptotic behavior for the ATM European call options in this mixed model was obtained in [381 using 
Fourier methods. We present, in Appendix [b] a probabilistic proof based on (4.1) and following an approach similar 
to that in T^. 

Proposition 4.1. In the setting of Section^ the at-the-nioney European call option price has the following asymptotic 
behavior: 



limt-^/^E{Sf - So)+ = SQaE*{W^)+. 



(4.2) 



Next, we give the second-order correction term for the at-the-money European call option price. As before, we 
change the probabilit y mea sure P* to P so that Xt = Z t + i7 + aW^ , with (^t)t>o a symmetric strictly F-stable Levy 
process under P (see (2.11)) and 7 defined as in (2.10). Recall also that, under both P* and P, W* is still a standard 



Brownian motion. We will also make use of the decompositions (2.14|-(2.17). The proof of the following result is 
presented in Appendix [BI 



Theorem 4.2. In the setting of Section^ the at-the-money European call option price is such that: 



lim 1 2 

t->-0 



-(t-i-E(.. 



So)+-aE*{W^ 



Y{Y 



— E* 



*\1-Y\ 



(4.3) 



Remark 4.3. As well known, the (1 
(25.6) in inn): 



y)-centered moment of a standard normal distribution is given by (see, e.g. 
21-^ 



-r 1 



Thus, the second-order asymptotic behavior of the ATM call option price (1.1 ) in short-time takes the form 

1 



So 



EiSt ~ Sor = dit^ +d2r 



(t^o), 



with 



di = 



d2 = -^r 



Y 

'-2 



Ca 



l~Y 



Y{Y-iy 



(4.4) 



(4.5) 



./27r V^r 

Intuitively, the first-order term only synthesizes the information about the continuous volatility parameter cr, while 
the second-order term incorporates also the information on the tail index parameter Y and the overall jump-intensity 
parameter C . However, these two-terms do not refiect the relative intensities of negative or positive jumps (controlled 
by the parameters G and M). This fact suggests that it could be necessary to develop a third-order approximation as 
described below. 



Remark 4.4. The proof of Theorem 4.2 (See Appendix [B]) also provides higher order asymptotics of the ATM call 



option price under the generalized CGMY model. In fact, as mentioned in the proof, the second integral in the 
decomposition of Bt is asymptotically equivalent to (7/2) ^/t, while the l ast in tegral is clearly asymptotically equi vale nt 



to 



-(T^/4)-\/<. Then, it remains to analyze the first integral At in (B.4), which in the proof of Theorem 



decomposed into /i(i), l2{t) and Iy,{t). For Ii{t), it can be shown that the second order te rm of Ji2{t,y) is 0{tr ) 



4.2 
1 



IS 



while the first order of Jx\{t,y) is 0{\/t). For /2(i), the first term in the decomposition (B.18) is o(-\/i), while the 



second order term is Oit^ ■^). Finally, for Iz{i), the second order of j'^(t,y) and jl^{t,y) is 0{t^ ^), while all the 
other terms in the decomposition of I^it) are of order o(-\/i). Therefore, as t — > 0, 



E{St - So)+ = + d2t^ + dst + o{t), 



1<Y<-^, 



So 

^ EiSt - So)+ - d,t"^ + d2t"-^ + daii-^ + 0(^5-^-), ^ < y < 2, 
bo 2 



(4.6) 
(4.7) 



where d^ can be explicitly derived. 
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The next proposition gives the small-time asymptotic behavior for the ATM Black-Scholes implied volatility, again 
denoted by a(t), under the generalized CGMY model. Unlike the pure-jump case, we can only derive the first order 
asymptotics using Theorem |4.2[ In fact, the first order term of the ATM call option price under the generalized CGMY 
model is the same as the one under the Black-Scholes model. The second order term of at requires higher order 
asymptotics of the ATM call option price. The proof is deferred to Appendix [B] 



Proposition 4.5. Under the exponential CGMY model (1.2) with non-zero Brownian component, the implied volatility 
a is such that: 

^W = -+^||^^r(i-|)t-¥+o(t-¥), t^o. (4.8) 



5 Numerical examples 

In this part, we assess the performance of the previous approximations through a detailed numerical analysis. 



5.1 The numerical methods 

Let us first select a suitable numerical method to compute the ATM option prices by considering two methods: Inverse 
Fourier Transform (IFT) and Monte Carlo (MC). 

Before introducing the IFT method, let us set some notations. The characteristic function corresponding to the 
Black-Scholes model with volatility S is given by 

2-t 



Vt m = exp I — [v +iv) 



The corresponding call option price at the log-moneyness k — \og{So/K) under the Black-Scholes model with volatility 
S is denoted by Cg^i^) ' that is, 



Let us also recall that the characteristic function under the generalized CGMY model is denoted by ipt (see (2.2)) 
and let us denote the corresponding call option price at log-moneyness k by C{k). The IFT method is based on the 
following inversion formula (see Section 11.1.3 in [5]): 



oo 



z^{k) C(fc) - C]^sik) = ^ y_ e-^^''Uv)dv, (5.1) 
where 

C.(^):-e- ^-(--)"57(--) . (5.2) 

iv(l + IV) 

In our case, we fix r = and, since we are only interested in ATM option prices, we set k = 0. In order to compute 
numerically the integral in (5.1), we use the Simpson's rule: 

p-i 



with A = Q/{P - 1), Vm = -Q/2 + mA, and wq = 1/2, W2e-i = 4/3, and W2i = 2/3, for £ = 1, ... , P/2. 
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We also consider a Monte Carlo method based on the risk-neutr a.1 opt ion price representation under the p robab ility 
measure P. Under this probability measure and using the notation (2.151 as well as the relations (2.14) and (2.17), we 
have: 



1) 



-Ut 



(1 



-Xt} 



the variables and 



Urp are 



which can be easily computed by Monte Carlo method using the fact that, under 
independent F-stable random variables with scale, skewness, and location parameters TC\ cos(7ry/2)|r(— F), 1, and 
respectively. Standard simulation methods are available to generate stable random variables. We consider the following 
set of parameters for the CGMY component: 

C = 0.5, G = 2, M = 3.6, Y = 1.5. 

Figure [T] compares the first- and second-order approximations as given in Remarks |3.3| and |4.3| to the prices based 
on the Inverse Fourier Transform (IFT-based price) and the Monte Carlo method (MC-based price) under both the 
pure-jump case and the generalized CGMY case with a = 0.4. For the MC-based price, we use 100, 000 simulations, 
while for the IFT-based method, we use P = 2^^ and Q = 800. As it can be seen, it is not easy to integrate numerically 



the characteristic function (5.2) since in this case T is quite small and, therefore, the characteristic functions ipT and 
if^ are quite flat. The Monte Carlo method turns out to be much more accurate and faster. 




Figure 1: Comparisons of ATM call option prices for two methods (Inverse Fourier Transform and Monte-Carlo 
method) with the first- and second-order approximations. MC-based price is based on 100, 000 simulations while the 
IFT-based method is based on the parameter values P — 2^^ and Q — 800. The parameter a in the generalized CGMY 
model is set to be 0.1. 



5.2 Results for different parameter settings 

Here, we investigate the performance of the approximations for different settings of parameters: 

1. Figure [2] compares the 1st- and 2nd-order approximations with the MC prices for different values of C, fixing the 
values of all the other parameters. In the pure-jump case, the 2nd order approximation is significantly better for 
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moderately small values of C, but for larger values of C, this is not the case unless T is extremely small. For a 
nonzero continuous component, the 1st order approximation is extremely bad as it only takes into account the 
parameter a. 

2. Figure [3] compares the 1st- and 2nd-order approximations with the MC prices for different values of Y, fixing 
the values of all the other parameters. In both cases, the 2nd order approximation is significantly better for 
values of Y around 1.5, which is consistent with the observation that |c?2| — >■ oo as F 1 or F 2. For a 
nonzero continuous component, the 1st order approximation is again extremely bad as compared to the 2nd order 
approximation . 

3. In the left panel of Figure |4j we analyze the effect of the relative intensities of negative jumps compared to 
positive jumps in the pure-jump CGMY case. That is, we fix the values M to be 4 and consider different values 
for G. As expected, since the first order approximation does not take into account this information, the 2nd-order 
approximation performs significantly better. 

4. In the right panel of Figure|4j we analyze the effect of the volatility of the continuous component in the generalized 
CGMY case. The 2nd order approximation is, in general, much better than the Ist-order approximation and, 
interestingly enough, the quality of the 2nd order approximations improves as the values of a increases. In fact, 
it seems that the 2nd-order approximation and the MC prices collapse to a steady curve as a increases. 




Figure 2: Comparisons of ATM call option prices with the short-time approximations for different values of the jump 
intensity parameter C. 
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Figure 3: Comparisons of ATM call option prices with the short-time approximations for different values of the 
tail-heaviness parameter Y . 




Figure 4: Comparisons of ATM call option prices with the short-time approximations for different values of G or M 
and different values of the volatility parameter a. 
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A Proofs of Section [3]: Pure-jump CGMY model 

For simplicity, throughout this section, we fix 5*0 = 1. 
Proof of Lemma 13.11 



From (1.6), we have 



/•CXD 

r^/^E{St - So)+ = t-i/^P* {Xt >E)^ t-^'^ / e-"P* {Xt > x)dx 

Jo 

e-*''"'"P*(i-i/^Xt > u)du. 



Next, using the change of probabihty measure ( |2.13| ), 



and, moreover, since a = 0, (2.11), and (2.14) 



>u-7ti-i/i'} 



) 



with Ut := M*Ul - G*U^ . By the self-similarity property (2.12), 



4Zi>ti-7ti-i/i'} 



Thus, plugging (A. 2) into (A.l) 



and, changing variables {y = u — ')t^^^/^^, the result follows. 



^ du, 



Proof of Theorem 1372 



To begin with, we assume that M — 1 = G, so that 7 = (see (2.10)) and, in light of Lemma 3.1 



-*''"-E (e~'" ■ ^-\{z,>,}] dv. 

The general case is resolved in Lemma [A. 1| below. Let 

D{t) := i-i/^E(5, - ^o)+ - E(^+), 

which can be written as 

D{t) := j e-*''"^E (e'*''"^i Mz,>v}) dv - E(Z+) 
+ (e-"* - 1)E(Z+) 



+ (e-"* - 1) ^ e-^'^'-E (e-'''"^^l{z,>,}) dt; - E(Z+) 
--■.Di{t)+D2{t)+D:i{t). 



We will show that 



t^/^-^Di{t) — >d + ri, as i ^ 0, 
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while it is clear that D3{t) = o{Di{t)) and t^/^^^D2{t) = o(l), as t 0. First, note that in light of ( |2.26 l 



Di{t) = E e 



'du ) - E(Z+) 



(^E(e-*^'"^^^) -E(e-*^'"(^^+^0)j -E(Z+) 
= (^e''*-E(e-*''"(^^+^0)^ -E(Z+). 



Thus, 



e'7* _ 1 /-"^ e-*'^''^ - 1 



+ / fi-i/Y ^^^t + Ui> v)dv - / ^-^P(Z+ + C/i < -v)dv 



Dii(t) 



Clearly, 

Next, for £)i3(i), note that 



-Dl2(t) 



Dii{t) rj, as i -J> 0. 



I5l3(t) 



¥{Z+ + Ui< -v) < P(C/i < -v) < 



E(£ 



and, since < — 1 < ye^ for y > 0, 



^i_i/y ■ + Ui< -v) < t^^^-'e^ve^'"'"-'^^ < e^ve-^l'', 

for t > small enough. The dominated convergence theorem then implies that 

-DisW ^ 0> as i ^ 0. 
Finally, let us analyze the term Di2(t). Changing variables [u — t^^^v), 



and, from (2.22) and (2.231, there exists, as shown in Appendix [C] a constant k < oo such that 

for any < t < 1 and u > 0. Hence, by the dominated convergence theorem, 

limDi2(t) = ^ {e-'' - 1) lim (^r^¥(^Z+ + Ui > t^'^^^uj) du. 
To find linit^o i"^P(Z+ + Ui> t-^/^u), note that 



-f2 
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Then, 



h ^ p((Af* + !)[/+ - (G* - l)C7f > t-^/^u, + t/f > o) 



M* + 1 



\ ^ ~ M* + 1 ^ - M* +G 

A,i(t)+/i,2(i). 



+ P ( f/+ > -U^ > — — 



By the independence of Ui and J/j^ and the estimate (2.20) 

/i,2W<P C/i+> 



M* + G* 



~ur> 



M* + G* 



(A.14) 



as t ^ Q. Note also that 



t I 1 - M* + 1 



1- 



M* + 1 



t Y u 
M*+G* 



and so, recalhng that y) denotes the density of —J/j^ , then using (2.20), Lemma 2.1 and the dominated convergence 
theorem. 



Urn = hm I I p{l,y)F { U+ > 



t-^'^u-{G* - l)y 



M* + 1 



li'— M'+G' I 



p(l,v) lim -P [/+ > '— 1, 

t^ot \ ^ - M* + l ) {y<h^} 

1 iC(M* + l)^ 



dy 



p{\,y) hm 



t-^o t Y{u - t^/^{G* - l)yY {y<w^} 
C{M* + l)^ 



(A.15) 



Similarly, 



/2.l(i)+/2,2(t). 



-?7f > > 



G* 



-c/r > tj+ > 

1 ^ - M* +G* 



Again, as in (A.14), 



'U[- > U+ > 



M* + G* 



(A.16) 



and since 



t Yu — M*y 
G* 



{y<'M^} - t 



G* 



{y<li4G^^} - t 



t Y u 

M*+G* 
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by (2.201, Lemma 2.1 and the dominated convergence theorem, 



, 1 ,^ , ^ f , f - t-'^/^u- M*y\ 

In. . lm.t-y^Ml,.)P (-C/f > ) \y^^^}dy 

, - t-^/^u- M*y\ 
p(l,,)ln.i-P (-U^ > 



C{G*) 
Yv} 



(A.17) 



Combining (A.8|, (A.9|, (A.ll), and ( A.14)-(A.17) imphes that 



limii/^^iZ?i(t) = -^((M* + 1)^ + (G*)^) (1 - e-") y-'^dv + r]. 
Finally, we use the following identity (see p. 84 in 'ST]): 

(e"" - i)v-^dv = r(i - y) = -yr(-y). 



This concludes the proof. 



□ 



Lemma A.l. If ^ ^ in (3.1), then 



1 1 



7 



(A.18) 

Proof: Without loss of generality, fix 5*0 = 1 and also assume that 7 > 0, the case 7 < being similar. Then, using 



(3.1) 



limi^/'*'-! {t-^'^¥.{St - 5*0)+ - E(Z+; 



limti/5--i (e-(7+';)t / ^ 

t-i-O 



+ limii/^-ie-(^+'')* / e 

t-fO 



"El 



7*1-1/'*' 



4Zi>d} 



:=i?ii(t)+I?i2(t). 



As in the proof of ( A.6 1, it can be shown that 



\unDii{t) =d + ri. 



t-i-O 



(A.19) 



For Di2{t), changing variables {u — t^/^ ^u), we have 



du — e 



9t {u)du, 



where 



Since Y £ (1, 2), it is easy to see that, for —7 < u < and < t < 1, 

\9i'\u)\ < P(Zi > -7, [7i > 0), \g['\u)\ < em{e'^^l^z,>-^}) 



{Zi>ti-i/^u}J-{c/i<o} 
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and 



'-{Zi>ti-i/i'«}l{£fi>o} ^ l{Zi>«}i 



L{Zi>ti-i/v-„}i{f;^<0} 



< e 



L{Zi>u}- 



By the dominated convergence theorem, it follows that 



\imDi2{t) 



(f{Zi > 0, [/i > O) + P(Zi > 0, f/i < 0)) du 



Combining the previous limit with (A.19I leads to (A.18I. 



(A.20) 

□ 



Proof of Proposition 3.5 



The small-time asymptotic behavior of the ATM call-option price CBs(t,cr) at maturity t under the Black-Scholes 
model with volatility cr and zero interest rates is known to be such that (fixing for simplicity 5*0 = 1) 



V27r 24v27r 



(A.21) 



see, e.g., [181 Corollary 3.4]). To derive the small-time asymptotics for the implied volatility, we need an analogous 
result to (A.21 ) when a is replaced by (T{t). To show such a formula, the following representation due [331 Lemma 3.1] 
will be useful 

CBs{t,<j)^F{<jVi) with F{e):= <S>'(-)dv=— exp{-vy8)dv, 

Jo v27r Jo 

together with the following Taylor expansion for _F at 6* = (see [33l Lemma 5.1]) 

1 „ 1 



m = 



/2^ 24^/2^ 

Then, since ^{t) — > as t — ^ (see, e.g., [311 Proposition 5]), we conclude that 



CBsit,CT{t)) 



^jt) ^1/2 ^{tf ^3/2 
V2^ 



24V2 



as i — ^ 0. 



(A.22) 



Returning to the proof of Proposition 3.5 by equating (3.4 1 and (A.22 1 and comparing the first order terms, 

V 27r 

and, therefore, 



<T(i) - %/2^E(Z+)t7"5 crii7"5^ t^O. 



(A.23) 



Next, set a{t) = a{t) — ait^ 2. By comparing the first and second order terms in (3.4) with the first term in (A.22) 
(noting that the second order term in (A.22) is o(t)). 



CT(-r) 



{{M - if - - (G + 1)^ + G^)t^ 



Hence, (j{t) — )■ as t 0, and moreover 

m - ((M - 1)^ - - (G + 1)^ + G^) t^'\ 



t 0. 



t -> 0. 



Combining (A.23) and (A.24) finishes the proof. 



(A.24) 

□ 
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B Proofs of Section [4]: CGMY model with Brownian component 



Proof of Proposition |4.1 



From (1.6), note that 



/■OO /"OO 

t-^/^E{St - So)+ = t-^/'^So¥*{Xt >E)= r^/^So / e-^P*(Xf > x)dx = Sq e-'^*''¥*{t-^/'^Xt > u)di 

Jo Jo 

Now for any u > and < i < 1, 



'{i-^l'^Xt >u) < P*(i-i/VVF; > u/2) +V*{t-^/^L* > u/2) 
:^F*{aWl > u/2)+¥*{t-^/^L; > u/2). 



(B.l) 



Clearly the first term in (B.l ) is integrable on [0, oo). To estimate the second term, applying the change of probability 
measure (2.131 and using the self-similarity property (2.12) of Z, we obtain 



'{t-'/^L; > u/2) = E (e-^'l^,.,,.^^>„/2-,yt}) = E ( 



{Zi>t 2 ~ y u/2-7r" Y } 



Pick 1 < q <Y and q' > 1 such that q ^ + q' ^ = 1, then by Holder's inequality and (2.26), 



E 



V {Zi 



>t2 Y} 



< E e 



-q'Ut 



1/9' 



1 1^ 

t2 YU 1 1 

Zi > — 7: jt 



1/9 



-It 



(E(e-«'^* 



1/9' 



YU 1 

1-2 ' 



1/9 



, /Y-l. 



'/ til-u 1 1 \ 



1/9 



< exp 



i 1_ 

t2 Y U 1 1 

Zl > ^t^-Y 



1/9 



since as given in (2.16), 77 > 0. If 7 < 0, then by (2.23), for any u > and < i < 1, 

1/9 



Zl > 



t 2 Y U 



1/9 



~ 1 1 



< 



8k 



' t2 Y u _ ~,j.l-J, 

^ 2 



It' 



< 



(16k)1/'3 

7/i'/9 



and thus. 



'/ t2 YU _ 1 



1/9 



< min 



(16k)1/« 



which is integrable on [0, +00). On the other hand, if 7 > 0, 



Zl > 



1 1^ 

i 2 Yu 



Zl > 



1 L 

i 2 1- u 



< 



16k 2^7 



7i'"^)l{«>2,Vi}+IP( 



^1 > 



1 1^ 

1 2 li 



i<27Vi} 



and thus 



1 i_ 

t2 Y U 1 

Zl > — T-i^"^ 



1/9 



< min 



/ (16K + 2^7^)i/n 
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which is also integrable on [0,+oo). Therefore, by the dominated convergence theorem, 

nOO 

\imt-^/^E{St - 5o)+ = Sn V*icrW^ > u)du = Soa¥.*{Wl)+. 

The proof is now complete. □ 
Proof of Theorem liTH 

For simplicity, we fix = 1. Recalling that Xt ^ + aW^ under P* and using (1.6), the self-similarity of W* , and 
the change of variable u — i~^/^a;, it follows that 

POO I-OO 

Bt:=t-^/'^E{St-So)+-aE*{W^)+= J e-^'^'V* (^aW^ > u - t'Htj du - J V*{aW^>u)du. 

Next, by changing the probability measure to P and using that L'^ = Zt + jt, Ut = Ut + r]t, and the change of variable 
y — u — in the first integral above, we get 



OO 



^ g-('?+7)t 



+ J (^e-'"-'^'y - Ij P*(crW^* > y)dy. 

The last term above is clearly 0{y/i) as i — 0, while the secon d term above can be shown to be asymptoti cally 
equivalent to (7/2)-\/t by arguments analogous to those of (A. 20). Thus, we only need to analyze the term in (B.2) 
that we denote At and that can be written as follows in light of the self-similarity property of Z and U : 

= --'"^ (^"^'^ r ^'^^(\.w,>,-*f -iza - U^w,>-y}yy) , (B.3) 

where we had denoted ry := r/ -|- 7. To analyze the asymptotic behavior of At, we decompose it into the following three 
terms: 



{w*>o.,crw^+tY-'2 Zi>a} 



aWT 



I {n<<yWl<-tY-2Z^} ^1 

(1 - 1 
e-*^^il^ 11,/ e-^'^dy 

■.^ hit)- hit) + hit). (B.4) 
We analyze each of three terms above in the following three steps: 
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Step 1. We first analyze the behavior oi Ii{t). Since {Zt)t>Q and (Wt*)t>o are independent, 
/i(t) = e-''*E 1 1 1 ^ e^VtaW^ 

' \ {W^>0,aW*+tY-2Zi>0} 



{Zi>-t2-Yy} 



dy 



.-fit 



JAt,y)e-^'y dy. 



Using that the distribution of Zi is symmetric under P (hence, EZi = 0), Ji{t, y) is then decomposed as: 



(B.5) 



Ji{t,y) = E 1, 11, 

:= Jii{t,y) + Ji2{t,y). 



-tYUl 



-tY(Ut+Zt) 



Vt 



tY-2Zl +tY-2E{Zll, 1 1 , 

' ' ' {Zl>t2-Yy} 



Let us first consider Ji2(t,y)- By (2.241, it follows that, for any < t < 1 and y > 0, 

t^~^Ji2{t,y) < Ayi-^, 

for some A < oo (see Appendix [C] for the verification of this claim). Moreover, for any fixed y > 0, 

/•oo poo 

t^~^Ji2it,y)=t^'^Y-2 I ^ upziu)du ^ t^^Y-2 / wpz{t^^^w)dw. 

Jt2^Yy Jy 



(B.6) 
(B.7) 



Using ( |2.25| ), 

for t small enough and all w >y. Therefore, by the dominated convergence theorem and, in light of (2.25), we get: 



t->o 



limt^-ie"''* / Ji2(t,y)e-^«-^-^ 



%/27rcr2 



dy 



(^limi2-ij^2(i,2/) 



dy 



oo / /"OO 



w^lirni2 y ipz{t2 Yw))dw 



V2^ 



fdy 



For Jii{t,y), note that 



oo / poo 



Cw ^ dw 



\Jy 



Jn{t,y) = t^-^E f l{z,>o} y^'''''^' (e-*^" - l) 



6 2^, C 



V2 



Y - 1 



.i-v 



e 2o 



%/2 



- 1 -t^-^E 1 1 

I {-t2 y'y<Zi<0} /,-. 



Ui+Zi 



(e-*""-l) 



(B.8) 



du 



tY 2 



\ 1_ 

tY 2 



1 1 
+ tY-2 



^1 > 0,[/i < u < Ui + Zi)du 
- l\Vi-t'^--^y < Zi <0,Ui + Zi <u < Ui)du 



l) (P(^i >0,Ui<u<Ui + Zi) 



t^^yy < Zi <0,Ui + Zi <u < Ui)]du 



l) {P(Zi >0,Ui<u<Ui + Zi) - V{-t^^-Yy <Zi<0,Ui + Zi<u< Ui)]dx. 
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Next, change variable back to x — t^^^u: 

^ /"OO 



-t's-Yy < Zi <0,Ui + Zi < r^x < Ui) ]du 



. {e-'^ -l){V{Zi>0,Ui<t-^x<Ui + Zi) 
-V(-t^-^y < Zi<Q,Ui+ Zi< r^x < Ui)^dx. (B.9) 

For t > and y > 0, set 

Ti{t,x,y) :=P(Zi > 0,Ui < t'^x < Ui +Zi), T2{t,x,y) :=P(-t5~^y < Zi < 0,Ui + Zi < t'^x < Ui). 



By (A. 12), there exists k> such that for any x > and <t <1, 

Ti{t,x,y) < V{r^x <Ui + Zi) < ktx^^. 



(B.IO) 



Hence, 



< e-'^H^-^ 



OO /'OO 



y-3 
< t~ 



Y-1 



Jo 

OO poo 



Ti{t,x,y)dx 



-dy 



e 2o 



(1-e '^)Ti{t,x,y)dx — dy 
^0 v27r(T^ 



OO /"OO 



e 2a 



(1-e '^)x dx =dy — > 0, 
Jo Jo v^na^ 



(B.ll) 



as i — ?> 0, since Y > 1. Similarly, using Ui — M*U^ — G*Ui and Lemma 2.1 for any < t < 1 and x,y > 0, we have 

T2{t,x,y) <¥{t-^x < Ui) 



2M* 



2G* 



t Y X 



< 2P f/+ > , 

\ 2{M* +G*) 

< 8Kt{M* +G*fx-^ . 



Therefore, 



rjt foo poo 



lim — 



2 Jo Jo 



T2{t,x,y)dx 



e " 



V27rcr2 



-dy = 0. 



(B.12) 



(B.13) 



For X < 0, since t/+ and -C/f are identically distributed, we proceed in the proof as follows: 

Ti{t,x,y) < P(c7i < t-i-x) = P{M*U+ - G*C7j- < r^x) 

<V{U+ < {2M*)-h-^x) +P(-i7f < {2G*)-h-^x) 

< 2P(C7+ < {2M* +2G*yh-^x) 
/ i_ \ 

t Y X \ 



< 2E(e-^i')exp 



2(M* +G*) 



(B.14) 
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which is again independent of y and integrable on [0,00) when multiphed by e ^ — 1. Moreover, 



dx 



0<t~ / {e-'' -l)Ti{t,x,y)dx 



" / I t-^x \_ ( t-^x 



r7+, t^{M* +G*f 



l~2tY{M* +G*) 
Hence, by the dominated convergence theorem, 



< e-^^T-i 



e-"^ - 1 



J -00 



Vi 



Ti{t,x,y)dx 



g-\/tyg-^ 



%/27rcr2 



-dy 



< 2E 



l-2t7{M* +G*) Jo V2 



(B.15) 



as t ^ since 2/r > 1 > (3 - r)/2, for 1 < F < 2. Similarly, since Ui + Zi ^ MU{ - Gt/f , it follows from ( |B.14[ ) 
that 



72(t, a;, y) < P(AfC/+ - GL/f < t^^x) < 2E(e~'^i' ) exp 



t Y X 

2{M + G) 



Therefore, 



lime-''*t^-i / / — 'L^T2{t,x,y)dx'—^^J^dy^Q. 



Combining ( |B.8| , ( |B.11| , ( |B.13[ ), ( |B.15[ ) and ( |B.16[ ), we obtain 

limi^~^/i(t) 
t— >o 



2/ 



Y-lJo V2na^ 



%/27rcr2 



(B.16) 



(B.17) 



Step 2. Next, we analyze the asymptotic behavior of l2{t)- Using the independence of {Zt)t>o and {Wf)t>o, 

1 - e^^""^!* \ 



./n V {Zi 



{0<aW^<~tY-2Zi} 

1 - 



E 



1)1 



1 - e-^^y ert, 



e 2' 



-:dy 



1 1 \ 1 - 



Jo ^ ~ ^ Vi 



By (2.23) and the symmetry of Zi, 



Zi<-t^-^y) 7= — =v{Zi>t^-^y) j= — <8Kt^^^y^^^ < Sny^^^ 



Vi 



Vi 



(B.18) 
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which, when muhiphed by exp(— y^/(2(T^)), becomes integrable on [0,00). Hence, by (2.241 and the dominated con- 
vergence theorem, 



hmi2 

t-^Q 



Zi < -t^ yy] ^ dy 



hm t 2 



V27rcr2 ■ 



hm t 



^i>t^ i-y)y^j^f=dy 



e 2^ 



\/2^ 



hmt2 ^ — (^2 yy) ^ y—- 



_ C 

To find the asymptotic behavior of the 
E((e 



y 

y n: — ^2 y 



first integral in (B.18), we decompose it 



{Zi<-i2-yy} 



as: 



: E( 1 



11-,/ e ^ 



(B.19) 



-E l, i_ 1 ~ , / 

J2i(t,y) + J22(i,2;). (B.20) 
For J2i(i,y), note that for any < i < {M* + G'*)~'*' and y > 0, 

J —00 

< f e~^P (a/* C7i+ - G* C7f < r ^ a;] dx 

J —00 

<2/_\-p(t7.-<^)o,. 

^ 'l-t^{M*+G*) 

which is independent of y. Since 1 — Y/2 < 1/2 < 1/Y, for 1 < F < 2, by the dominated convergence theorem, 



< t^-i2E(e-^i ) \ - 

^ ^ 1 - tT {A-I* + G* 



_ y 
e 2^ 

^y27ra 



.dy — > 0. 



(B.21) 



We 



further decompose the second term J22{t,y) in (B.20) as: 

-t^E(Uil^ 1 1 _ , 

V {Zi<-t2-Yy,Ui>0} 

■.^jg\t,y)-jg\t,y). 



(B.22) 
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Since 1 - r/2 < 1/2 < l/F, for 1 < F < 2, it is easy to see that 



Moreover, 



(B.23) 



tYUl 



{Zl<-t2-Yy^Ul>0} 



</ (1 -e""")?!^;/! > t-^wjdu;. 



Hence by (B.12| and the dominated convergence theorem, 
Jo vt v27ro"^ 



°° 1 - e-™ 



w 



dw / y^=,dy ^ 0, (B.24) 



"72 



as t ^- 0. Combining (B.19I, (B.21I, (B.23) and (B.24), we obtain 



Y C 

Umi¥"^/2(i) = — , 



(B.25) 



Step 3. We finally analyze the behavior of Isit). Note that 



' ' \ {Zl>t2-Yy} 



1 



1 _ev^ye-*^^i 
1 - e^^-y e"^ 



72^2 



{Zi>t2-V-y} ^ 



I e^^ dy 



We first investigate the asymptotic of J3i(t, y) by decomposing it as 

J3i(t,2/)=Ef(e-**^^-l)l 



V27r(72 



(B.26) 



(B.27) 



By (2^1, it is easy to see that J^f{t,y) < Snt^^^y^^, for any < t < 1 and y > 0. Hence, by {'L24] and the 
dominated convergence theorem, 



lim — 



vt (2).. A~e^*y 



2 Jo 



Vi %/27rcr2 



dy 



V27rcr2 



dy. 



(B.28) 



Also, ^3i''(i,y) can be further decomposed as: 



E 



-tYUl 



1)1 



{Zi>th-^y,Ul>0} 



e"P Zi > t^-Yy^Ui < -rYu]du 
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For any u > 0, y > and t > 0, let T^it, u, y) :— f(^Zi > ^ y,Ui < —t y u^. It is easily seen that 

i_ 

m,u,y) < p(c7+ < +p(-c7r 



< 



< 2Ee^i" exp 



-t YU\ 

2M* J 

2{M* +G*) 

1^ 

— t Y U 



— t Y u 



2G* 



2{M* + G* 



Moreover, 



< E (1 



-tYUl 



)1 



{Zi>t5 ^y.Ui>0} 



< tYE\Ui\ 



(B.29) 



(B.30) 



Hence, for any y > and < < < 2-^ (M* + G*)"^, and since 1 - y/2 < 1/2 < l/F, for 1 < F < 2 



<i^"'2Ee^i" / exp 







-t Y u 



2{M* +G*) 



4E e^i k^(M* +G*) 
-1 — V — y ht^ + ^" 



1 - 2tY{M* +G* 



^EU7i — ^ 0, as i ^ 0. 



(B.31) 



Since both control functions in (B.29) and (B.30) are independent of y, combining (B.28) and (B.31), and by the 
dominated convergence theorem, 



(B.32) 



Next, we decompose the quantity Ji2{t,y) defined in (B.26) as: 



J32{t,y) = E 1, 1 1 

•^^^ ' \ {Zi>t2-Yy} 



tY-2Zl 



V {Zl>t2-Yy}J 



(B.33) 



similar to (B.8 ) 



Note that ,1^2 V'jy) is the same as Ji2{t,y) in (B.6), and thus the corresponding integral has an asymptotic behavior 



limi^-^e-'** rift^-iZil . . \e^'y^^dy=^ 



V Ina-^ 



(B.34) 



Next, we further decompose J^{t,y) as: 



= t^"5E 1^ 11,, 



Zi+Ui 



l)du 



t^^h [ (e-*^" - l)p(Zi > ti-^y,Ui <u<Zi + uAdu 



Vi 



ViJo 



(e-^ - l)PUi > t2-Yy^Ui < t-Yx < Zi + UAdx 



- l)v(Zi > t2-Yy^Ui < f-Yx < Zi + Ui]dx 
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Note that for a; > 0, 



while for a:: < 0, 



' Zx > t^^^y.Ui < r^x <Zi + < v(t^^x < Zi + Ui 



Using the estimates (B.IO) and (B.14|, a proof as in getting (B.ll) and (B.15) gives 



Combining (B.32 1, (B.34 1 and (B.35 1, we have 



C 



e 2<7 



.dy = 0. 



t^o Y{Y-l)Jo 



y 



1~Y 



V2 



(B.35) 



(B.36) 



TTCT^ 



Finally, from ( |B.4[ ), ( |B.17[ ), ( |B.25[ ) and ( |B.36[ ), and since 1 - y/2 < 1/2, for 1 < y < 2, we obtain (|43|, therefore 
finishing the proof. □ 



Proof of Proposition |4.5[ 

When the diffusion component exists, [551 Proposition 5] implies that d{t) -> cr as i — > 0. In particular, att^^"^ — ?> as 
i — > and, thus, (A. 22 1 above still holds. Let a{t) := at — cr, then a{t) — > as t — > 0, and ( |A.22 ) can be written as 



CBs{t, ait)) = + Mil/2 _ ^(^^3/2 ^ o / (a(i)t^/')') - + ^t'/^ + 0(i3/2). (g 



By comparing (4.4)- (4.5) and (B.37I, we have 



y(y-i)0FH ^ 



V2^ 



and, therefore, 



The proof is now complete. 



a{t) 



y(r-i) H 2 



t ^ 0, 



t ^ 0. 



□ 



C Additional proofs 



Verification of (2.8 1 



By the very definition of P* and (2.2 1, we have 



E*(e™^') = E(e(™+1)-^*) = ~ i) ^ gtCr(-Y)((Af-i(n-j))^+(G+»(u-i))^-M^-G^)+»c(«-»)t-4(n-»)2t 

_ gt Cr(-y)((M*-iM))^ + (G*+m)"^-M^-G'')+i(c+(j2)tit+ct-2^t+i^t 



with M* = M - 1 and G* = G + 1. Next, using ( |2.3| ), we clearly have 

<Gr(-r)(^M^ - G^) + ct+ta^/2 = -tCT{-Y) {{M - 1)^ - (G + 1)^) , 



and thus. 



_ Cr{~Y)[{M' -iu)^ + {G' +iu)^ -M'^ -G'^)+ic' ut- !!^t 
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with c* = c + cr^ . 



□ 



Verification of (2.10). 



Using (2.91, 



EXi = b+ / xi>{dx) = b*+ / x{i) - i^*){dx) + / xi>{dx) 

J{\x\>l} J\x\<l J\x\>l 

= c*+ f x{i>-iy*){dx)-CYT{-Y){{M*f-^ -{G*f-'). 



On the other hand, 



x(e^^*" - 1) dx+ / x(e-'^*" 

^ J-oo 



1) , ., ax 



oo 



E 



C{M*) 



/ a;"-^e'^^*-dx - V ^^^^ / 



n—l m — 1 

= -C{M* f-^r{l -Y)+ C{G* f-^T{l - Y). 



Hence, EXi = c* and (2.10) foUows. 



Proof of Lemma 12. It 



By (2.22), for any < i < 1 and w > 0, we have 



> t 



< 



1 /2C 



> t v] ('l{t-i/i'„>^} + l{t-i/i'i,<Ar}) 



t V 1^^^ l{i-t-i/^>Ar} + JpTjY:^Mt-^/Yv<N} 

< {2CY-' + N^) v-"" . 
The result follows by setting k = 2CY-^ + . 



Verification of (A. 12). 

Note that 



UI* 



+ r 



-Ui> 



AG* 



Using Lemma [O] there exists a constant < k < oo such that 



4Af* 



+ r 



46* 



< 2K(4Ar +AG*Y u 



for all w > and < i < 1. Clearly, (2.23) implies that, for any < t < 1 and w > 0, 



-P > ^ ] < 22^+1^^-^ < 32kw-^, 
t \ 2tY - 



□ 



□ 



and (A. 12) follows. 



□ 
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Verification of (B.7| 



In light of (2.24), there exist R> and H > 0, such that for any u> 

pz{u) < Ru-''-\ (C.l) 



Now for Ji2{t,y), using (C.l), for any < i < 1 and y > 0, 



^Ji2it,y)=t^ 2 I ^ upz{u)du 

Jt2 Y y 



i(y-2)(i-^)l 11 / upz(u)du 



t2 Yy 

^(t'-5)(i-^)i 11 / unz(u)du+ / upz(u)du 

OO 



<t(T-m-y)i 11 / Ru-'-du 

+ i(^-^)(i-ni^^i _ 1 ^^^^ (^J^ Ru~^du + HnZi > t^-^y)) 



OO / TT \ Y-1 



and (B.7) follows. □ 
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